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ABSTRACT 

Orthogonal Frequency Division Multiplexing 
(OFDM) has been recently applied wide in wireless 
communication systems, as a result of its high rate, 
transmission capability with high information 
measure, potency and its strength to multipath delay 
.Channel estimation is an important drawback in 
OFDM system. Channel estimation forms the center 
of any orthogonal frequency division multiplexing 
(OFDM) based mostly wireless communication 
receiver. Frequency domain pilot aided channel 
estimation techniques are either least squares (LS) 
based or minimum mean square error (MMSE) based. 
LS based techniques are computationally less 
complex. Unlike MMSE ones, they do not require a 
priori knowledge of channel statistics. However, the 
mean square error (MSE) performance of the channel 
estimator incorporating MMSE based techniques is 
better compared to that obtained with the 
incorporation of LS based techniques. We propose an 
alternative optimization algorithm which uses 
frequency-domain block-type training symbols, and 
prove that this algorithm always converges, at least to 
a local optimum point. Then, we calculate the 
Crame 'r - Rao lower bound, and show that the 
proposed estimator attains the bound. 
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I. INTRODUCTION 

Orthogonal frequency division multiplexing (OFDM) 
is widely recognized as one of the key techniques for 
high data rate communications in wireless networks. 
Due to its capability to mitigate multipath effect and 
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eliminate inter-symbol interference, OFDM has been 
utilized for current technologies. Despite its 
impressive qualities, OFDM signals suffer from large 
envelope fluctuations producing a high Peak-to- 
Average Power Ratio (PAPR). This leads to an 
unwelcome tradeoff between linearity of the 
transmitted signal and the cost of the High Power 
Amplifier (HPA). A number of approaches have been 
proposed to cope with the SISO-OFDM PAPR 
problem. First, techniques based on the channel 
coding transmit only the code words with low PAPR. 
Such coding techniques offer good PAPR reduction 
and coding gain. The critical problem of coding 
approach is that for an OFDM system with large 
number of subcarriers, either it encounters design 
difficulties or the coding rate becomes prohibitively 
low. Phase rotation is another approach to reduce 
PAPR, including selective mapping (SLM), and 
partial transmit sequence (PTS). It generates a set of 
sufficiently different candidate data blocks, all 
representing the same information as the original data 
block, and selects the most favorable block for 
transmission. 

Amplitude clipping may be one of the most effective 
solutions when the number of subcarriers is large. It 
clearly removes the amplitude peak, and does not 
introduce redundancy and power increase. Clipping 
the high PAPR signal by a limiter can be a very 
promising approach to reduce the cost of having 
highly efficient HPAs. Furthermore, clipping the high 
PAPR signal increases the battery life, thanks to the 
resulting increased power efficiency of the HPA. 
Estimating the nonlinear parameter at the receiver 
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side is important, especially in IoT-based MTC 
because there are a lot of sensors that need to send 
some information occasionally. Then, at each time of 
transmission, the receiver needs to estimate the CA to 
have an updated value which will be deployed in the 
detection stage. Those CAs for different nodes may 
change during time because of environmental 
conditions such as temperature, humidity, and 
unstable power supply. Furthermore, for low cost 
devices, the manufacturing variations from unit to unit 
may be fairly large and it would be too expensive to 
measure and calibrate each unit separately. Moreover, 
electronic devices and components suffer from 
process of aging, so even if this nonlinearity 
parameter was known at the beginning, it would not 
later. Therefore, these nonlinearity parameters will 
become unknown to the transmitter during time and 
hence there is an essential need to estimate them at the 


receiver side occasionally. In this paper the CA is 
assumed to be a priory unknown to the transmitter 
(and the receiver), and the receiver uses an alternating 
optimization algorithm to give the jointly optimal 
estimate of the channel taps and the CA. Once the CA 
and the channel have been estimated, the system uses 
them to detect the transmitted symbols by deploying 
the iterative detection method proposed in. 

OFDM reduces the effect of the Multi-Path fading by 
converting a frequency selective channel into a 
Parallel collection of frequency flat sub-channels. We 
investigate channel estimation technology in OFDM 
system. We mainly investigate the channel estimation 
technology in OFDM system. 



Fig: 1 System for channel estimation 


II. SYSTEM MODEL 

Consider the Fig: 1 an OFDM system with N number 
of subcarriers, in which S= [s 0 , s x , s 2 ,... ■ %-i] r is the 
frequency domain symbol vector selected from a 
constellation such as QAM. The time-domain symbols 
area unit obtained by taking the Inverse distinct 
Fourier remodel (IDFT) from the frequency-domain 
symbols as: 

1 -j2nn 

*n=-jjj 2fc=0 s k e ^~ , n = 0, 1., N-l (1) 


Fir) = { 


r, r < A 
A, > A 


0(r) = 0 


(3) 

(4) 


Where r is the magnitude of the limiter input signal, 
and A is the CA. Combining AM/AM and AM/PM 
conversion characteristics, g(u; A) as a function of 
complex scalar u (and parameterized by A) can be 
written as, 


giu-,A) = | 


u, \u\< A 
A e i ar g ( u ) ; \ u \> A 


(5) 


Rewrite this equation in matrix for as: 

x = F"s (2) 

Where ‘F’ is the Nx N unitary discrete Fourier 
transform matrix. The channel is slow fading with L+l 
taps (L«N), denoted as h = [h 0 , h lt h 2 ,.... h L ] T , g(. ; 
A) is the limiter non linearity with the following 
amplitude modulation and phase modulation 
conversion characteristics: 


The output of the limiter is z = g(x; A), in which g(: 
;A) is taken element-wise. However, it is difficult to 
directly work with the output of the limiter, where the 
output of the limiter can be represented in a linear 
fashion by introducing N augmented binary variables 
c n for n G {0,1,... .,N — 1} indicating whether the 
sample at time n has been clipped (c n = 1) or not ( c n 
= 0), i.e. 


c n 


1, r n > A 
0, r n < A 


( 6 ) 
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Which lead to , 

z n=(l- c n )x n +Ac n e )0n n— 0,1,...,N-1, 

(?) 

It can be represented in vector form 

z = (1 — c) O x + Ac O e- 70 (8) 

To remove the Inter-Symbol Interference (ISI), a 
Cyclic Prefix (CP) with a length Lcp (> L) is pre¬ 
added to the time-domain symbols at the transmitter 
and is removed at the receiver. After the process of 
adding and removing CP, the matrix form time-domain 
representation of the OFDM transmission can be 
written as: 

u =Hz + w (9) 

Where H is an N x N circulant matrix whose first 
column which represents the circular convolution 
operator, and ‘w’is a zero-mean circularly symmetric 
complex Gaussian noise vector, Taking the DFT from 
both sides of (9), we obtain the frequency-domain 
representation of the OFDM transmission as: 

y = Fu 

= y/N D H Fz + w (10) 


Where D H is a diagonal matrix with the N-point 
DFT of h as its diagonal elements, w is the DFT of 
time - domain noise vector. 

III. COMBINED CHANNEL & CLIPPING 
AMPLITUDE ESTIMATION 

Here propose an alternating optimization algorithm to 
estimate CA. To do this, we use frequency-domain 
block-type training symbols. Once we have the 
estimates, we are able use them to detect the 
transmitted symbols in the subsequent OFDM blocks. 

The intuition behind this algorithm is that CA is a 
slowly time-varying parameter (much slower than 
channel variations). Moreover, wireless channels are 
usually slowly time varying, so the block-fading 
channel model, which remains the same during the 
transmission of several OFDM blocks is reasonable. 
Using frequency-domain block-type training symbols, 
the problem of joint Maximum-Likelihood (ML) 
estimation of channel and CA can be formulated as the 
following Least- Squares (LS) using (10): 

m i n J 4>o,/? ||y — VND H Fz\\ 2 (11) 



Fig: 2 Basic idea for clipping amplitude estimation 


A. Clipping Amplitude Estimation Given the 
Channel 

First solve (11) for A given h. To do so, we first sort 
the elements in r, and construct a new vector denoted 
as f. Note that, r. = Ps r, in which Ps is the sorting 
permutation matrix. Using the sorted vectors, the 
estimate of CA at the i th iteration can be written as, 


O 

In which x = Psx, c = Psc, and cp = Pscp. Note 
that c is constant within each interval ?k— 1 < A < 
rk, k = 0...N -1. 



_ tttfrnW 

l-e~Y 2 )\H m [nW 


+ W 0 ’ 


0<n<N-l 


is the MMSE based estimator. 


A = arg min^ >0 ||y — y[N D l H FPJ z 
Where 


( 12 ) 


z = x 0(1— cj + AeJ® O c 

(13) 


To find the global minimizer, in each interval, we need 
just check the point at which the derivative is zero and 
the rightmost comer point. Therefore, we have 

Algorithml. Estimation of Clipping Amplitude (A) 
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1. Inputs: 

y, r and exp(j0) 

F and D H 

4 fork=0, 1, ....N-l 

2. Initialize: 

[r, P s ] = sort(r) 

e j0 = p s6 j0 

B = Vnd h fp s t 
f — 1 = 0 

3. for k =0 to N-l do 

A X — aH m[n] _ 

k (l-e-Y 2 )\H m [n\\ 2 +N 0 

5. A = A k 

B. Channel Estimation given Clipping Amplitude 

Here, we solve (11) for h given A. We can rewrite (11) as: 
minjy- VNdiag(Fz)Fh|| 2 (14) 

Therefore, the LS channel estimate can be computed as: 

hLS 

— (f h diag(Fz) H ) diag(Fz)F) f F H diag(Fz) H y 

(15) 

Where (.) + denotes the Moore-Penrose pseudo-inverse. 

C. Initialization of the Alternating Algorithm 

Two alternative initialization strategies to the 
alternating optimization algorithm as follows: 

1. Initializing by the Channel: Since the value of 
A is unknown at the beginning, channel can be 
estimated using the un-clipped version of 
transmitted time-domain symbols. It is equivalent 
to putting x instead of z in (15), and using s = Fx. 
Therefore, the initializing channel estimate is: 

C®= _ 

-!= (F H diag(s) H diag(s)F) 1 F H diag(s) H y 

(16) 

2. Initializing by the Clipping Amplitude: By 
substituting (15) into the cost function in (14), the 
resulting LS error for A is: 

a(A) = y H (I -T)y (17) 


Based on which initialization is used, we have two 
alternating algorithms given as Algorithms 2 and 3: 


Algorithm 2 Alternating Optimization with 

Initializing Channel 

1. Inputs: 

y, r and exp(j0) 

F 


4 for k =0, 1, ....N-l 
2. Initialize: 

[f, P s ] = sort(r) 

e J0 = Ps6 j0 

To) 


D^ = diag (Fh ° ) ,h ° is given by(21) 

r — 1 = 0 
i=l 

3. while(convergence criteria not met) do 

4. B = y/N D^FPj 

5. calculate A 1 " 1 ^ and using algorithm 1 

6. z = x © (1 - 4(;-i)) + / l_1) e j0 O 4(;-i) 

7. V = diag (PP/z)F 

8 . R h = E{VV h } 

(0 


9. hr 


— Rh 


10. D® = diag(Fh (i) ) 

11. i = i + 1 

12. end while 




4s (0 


Algorithm 3 Alternating Optimization with Initializing 

CA 


1. Inputs: 

y, r and exp(j0) 

F 

c' k fork=0, 1, ....N-l 

2. Initialize: 

[r, P s ] = sort(r) 

ei 0 = P s ei 0 

3. while(convergence criteria not met) do 

4. V = diag(FP s T z (i " 1) )F 


5. h 


(i-i) 


vrnmse 


— Rh 


^ + - 2 (=h) j 


h ls( 1 ~ 1 ) 


6. 1 ' ) = diag(Fl4 ^) 

7. B = Djj _1) FP s T 

8. calculate A® and k® using algorithm 1 

9. i = i + 1 

10. end while 
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Therefore, the objective function is decreasing at each 
iteration and eventually converges to a local minimum. 
The same argument is also valid for Algorithm 3.For 
Symbol Error Rate (SER) simulations use the iterative 
detection algorithm introduced in, but with the 
estimated channel and CA as follows: 

Algorithm 4 Iterative Detection Algorithm 

1. Inputs: 

y, F, and N q 
- (0 

hmmse and A calculated by Algorithm 2 (or 3) 

2. Initialize 

^H=Aia (Fh) 

k= \-e~ A2 + A Jj erfc(A) 

= 0 

3. for i=l to N q do 

4. s® = 0n 1 y ~ d a-1) )> 

5. x® = F h s® 

6. d ® =F(g(x® ;A) - x®) 

7. end for 

V. RESULTS 



LS & MMSE performance of the CA and 
channel estimation using Algorithms 2 and 3, 
when L + 1 = 7, 

N = 64 and CL = 1 dB. 


IV. CRAMER-RAO LOWER BOUND 

The Cram'er Rao Lower Bound (CRLB) expresses a 
lower bound on the achievable variance of unbiased 
estimators. The ML estimator asymptotically achieves 
the CRLB under some regularity conditions. The 
logarithm of the Probability Density Function (PDF) 
of the frequency domain observation vector y given in 
(10) can be written as: 

logp(y:0) =-j ||y — ViVZ) h Fz|| + constant 

Where constant comprises the terms which are 
independent of the estimation parameter vector 0 the 
log PDF is a differentiable function of ‘h’ in the whole 
parameter space. 

1 

CLRB(A ) =--— r - 

v J q-2P H C^P 

And 

1 

CLRB(Ji) = Cf 1 +--— z — CP^Ppn CP H 

W 1 q — 2P H C^ 1 P 1 1 



LS and MMSE performance of the chamr 
estimation using Algorithms 2 and 3, when L+l 
7 N = 128 and 
CL = 3dB. 
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SER performance of the iterative detectic 
(Algorithm 4) vs. SNR, when L + 1 = 7 N =6 
and CL = 3 dB. 


SER performance of the iterative detectic 
(Algorithm 4) vs. SNR, when L + 1 = 7 N = 12 
and CL = 1 dB 




SER performance of the iterative detection 
(Algorithm 4) vs. CL, when L + 1 = 7, N=128 
and SNR = 20 dB 

VI. CONCLUSION 

In this we, studied joint maximum-likelihood 
estimation of channel and clipping level at the receiver 
side in future IoT-based OFDM networks,. In 
particular, we have proposed two alternating 
optimization algorithms, And also computed the 
theoretical lower bounds (CRLB) on the performance 
of these estimators, and showed that they attain these 
lower bounds. Next, we have combined the channel 
and the CA estimates with the iterative detection 
method from to perform symbol detection at the 
receiver. Finally, we have showed by simulations that 
the performance of the iterative detection technique 
victimization the planned algorithm is sort of identical 
because the one amongst the case that the receiver 
have genie-aided data of the channel and CA. 


SER performance of the iterative detectic 

(Algorithm 4) vs. CL, when L + 1 =7, N = 51 

and SNR = 20 dB 
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